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. ^ 1 INTRODUCTION 



ABSTRACT 

We have made a series of snap-shot observations of 37 polars using XMM-Newton. 
We found that 16 of these systems were in a low, or much reduced, accretion state. 
Of those, 6 were not detected in X-rays. This suggests that in any survey of polars, 
around half will be in a low accretion state. We tested if there was a bias towards 
certain orbital periods: this is not the case. Of the 10 systems which were detected at 
low, but significant rates in X-rays, 8 showed significant variability in their X-ray light 
curves. This implies that non-uniform accretion still takes place during low accretion 
epochs. The bolometric luminosity of these systems is ~ 10'^'^ ergs s~^, two orders of 
magnitude less than for systems in a high accretion state. The X-ray spectra show no 
evidence of a distinct soft X-ray component. However, the X-ray and UV data imply 
that such a low temperature component exists: its temperature is low enough for its 
flux distribution to move outside the bandpass of the X-ray instruments. 

Key words: Stars: individual: ~ UZ For, CP Tuc, FH UMa, CV Hyi, AR UMa, V834 
Cen, UW Pic, V393 Pav, V4738 Sgr, RX J0153-59, MN Hya, ST LMi, RX J1313-32, 
V1033 Cen, RS Cae, MR Ser - Stars: binaries - Stars: cataclysmic variables - X-rays: 
stars 
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Polars or AM Her systems are accreting binary systems in 
which material transfers from a dwarf secondary star onto a 
magnetic {B ~10-200MG) white dwarf through Roche lobe 
overflow. They vary in brightness over the binary orbital 
period, which is synchronised with the spin period of the 
white dwarf. Unlike non-magnetic accreting binaries they do 
not have a disk and hence a reservoir of material. Therefore 
when mass transfer ceases, accretion stops. Many systems 
have been found to have episodes in which the accretion 
rate is much reduced, or accretion has ceased completely. 

The study of Lamb & Masters (1979) showed that when 
in a high accretion state {rn^l g s~^ cm~^), the density of 
the accreting gas in the accretion column is sufficiently high 
to maintain a colfisionaf timescafe shorter than the coof- 
ing timescafe and the flow dynamicaf timescafe. The flow 
is hydrodynamic, and an accretion shock forms near the 
surface of the accreting white dwarf. The shock heats the 
accreting materiaf to temperatures of severaf tens of keV 
and thereby emits keV X-rays (see Wu 2000 and references 
therein). There is afso a softer X-ray component, due to 
heating of the surrounding photosphere of the white dwarf 



by the hard X-ray irradiation or by bombardment of dense 
bfobs in the accreting materiaf. 

In contrast, in polars with fow accretion rates (m< 0.1 
g s^^ cm^^) collisions can be so inefficient that not only the 
efectrons and the ions become thermaffy decoupfed but afso 
ion-ion coffision cannot maintain the flow in hydrodynamic 
equifibrium. It was suggested that the kinetic energy of the 
accreting gas is refeased via Coufomb coffisions of charged 
particfes in the photosphere of the white dwarf (Kuijpers & 
Pringfe 1982; Thomson & Cawthorne 1987; Woefk & Beuer- 
mann 1992, 1993). The resufting spectrum woufd be signif- 
icantly different to that found in the high accretion state. 
This has been confirmed in X-ray observations of pofars seen 
in high and fow accretion states, see eg the fow-state spec- 
tra of ST LMi (Ramsay, Cropper & Mason 1995) and AM 
Her (de Martino et af 1998). In particular, the soft X-ray 
component seen in high accretion states is absent. 

However, even during epochs of fow accretion rates, 
the flow is very unsteady, with flares being seen on short 
timescafes and in various wavefengths. In the opticaf band, 
AM Her has exhibited flares in the fow state (eg Bonnet- 
Bidaud et af 2000). Ffares are afso seen in X-rays: in the 
pofar UZ For an accretion event was observed using XMM- 
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Newton at the phase when the primary accretion region was 
in view (Still & Mukai 2001, Pandel & Cordova 2002). The 
X-ray spectrum was consistent with a single temperature 
plasma of temperature ~4keV: no distinct soft X-ray com- 
ponent was evident. 

We have undertaken a survey of 37 polars using XMM- 
Newton (cf Ramsay & Cropper 2003 for a preliminary 
overview of this survey). For details of those systems ob- 
served in a high accretion state, see Ramsay & Cropper 
(2004) and references therein. Our survey also gives us an 
indication of how often polars go into a low/off accretion 
state. This paper reports on those polars which were in a 
low state, in particular their number, light curves and their 
spectra. 



2 OBSERVATIONS 

XMM-Newton was launched in Dec 1999 by the European 
Space Agency. It has the largest effective area of any X-ray 
satellite and also has a 30 cm optical/UV telescope (the Op- 
tical Monitor, OM: Mason et al 2001) allowing simultaneous 
X-ray and optical/UV coverage. The EPIC instruments con- 
tain imaging detectors covering the energy range 0.15-lOkeV 
with moderate spectra resolution. The OM data were taken 
in two UV filters (UVWl: 2400-3400 A, UVW2: 1800-2400 
A) and one optical band {V band). The observation log is 
shown in Table 1. The data obtained using the RGS instru- 
ment (den Herder et al 2001) were of very low signal-to-noise 
and are therefore not discussed further. 

The X-ray data were processed using the XMM-Ncwtoii 
Science Analysis Software (SAS) v5.3.3. For the EPIC de- 
tectors (Striider et al 2001, Turner et al 2001), data were 
extracted using an aperture of 40 centered on the source 
position. Background data were extracted from a source free 
region. The background data were scaled and subtracted 
from the source data. The OM data were analysed in a sim- 
ilar way using omichain and omfchain in SAS v5.4. Data 
were corrected for background subtraction and coincidence 
losses (Mason et al 2001). We show in Table 1 the mean V 
mag determined using the OM for each source at the time 
of our observation. 

For a typical polar in a high accretion state (eg V347 
Pav, Ramsay et al 2004), we estimate that the count rate 
will be a factor of ~4 greater using XMM-Newton compared 
to ROSAT. Every source shown in Table 1 was found to have 
a X-ray count rate which was significantly lower than when 
it was observed in the ROSAT all-sky survey with the excep- 
tion of AR UMa and ST LMi which wore not detected and 
MR Ser when it was in a low state. However, we extracted 
two pointed ROSAT observations of AR UMa and found no 
significant detection in either. Further, it has been found to 
show a high X-ray state and an optical range of ^14.5-16.5 
(Schmidt et al 1999) while ST LMi was found to be bright 
in X-rays using EXOSAT (Beuermann, Stella & Krautter 
1984). MR Ser was also detected using EXOSAT (Angelini, 
Osborne & Stella 1990). Using the highest reported count 
rate and an appropriate model, the expected count rate for 
MR Ser in a high accretion state using the XAIM-Ncwton 
EPIC pn detector is 4.2 ct/s. We conclude that all 16 sys- 
tems in Table 1 were in low or low/intermediate accretion 
states. 
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Figure 1. The phased EPIC MOS2 data (0.15-lO.OkeV) of V834 
Can phased on the ephemeris of Cropper, Menzies & Tapia (1986). 

Of these systems, 10 (UZ For, MN Hya, RX J1313-32, 
CP Tuc, V834 Cen, UW Pic, V393 Pav, ST LMi, V1033 Cen 
and MR Ser) were detected at low, but significant, levels. 
We discuss these systems individually in the next section. 
We do not discuss UZ For since these XMM-Newton data 
have been discussed in detail by Still & Mukai (2001) and 
Pandel & Cordova (2002). 



3 SYSTEMS SHOWING X-RAY VARIATION 
3.1 V834 Cen 

Sambruna et al (1994) presented multi-epoch observations 
of V834 Cen using EXOSAT and found that the shape of its 
low energy light curve (0.02-2keV) varied significantly from 
epoch to epoch. Our EPIC M0S2 X-ray fight curve (Figure 
1) shows significant flux at all orbital phases and with signif- 
icant modulation. (Both the EPIC MOSl and pn data were 
in timing mode and are therefore not considered further) . It 
is most similar to the 1986 EXOSAT light curve (Sambruna 
et al 1994). However, there is no evidence for an absorption 
dip at <j)=0.0, even at lower energies, although this maybe 
simply due to the relatively poor signal to noise (the ac- 
cumulated error on the ephemeris of Cropper, Menzies & 
Tapia 1986 is 0.08 cycles). 

We show in Figure 2 the X-ray spectrum taken using 
the EPIC MOS2 detector. We fitted the spectrum firstly 
using the stratified accretion column model of Cropper 
et al (1999). We fixed the mass of the white dwarf at 
0.66Mq (Ramsay 2000) and the specific accretion rate at 0.1 
g s~^ cm^'^ (appropriate for a system in a low accretion 
state). We find a fit with xS=l-49 (42 dof). There is no 
evidence for a distinct soft X-ray component. 

We then fitted the spectrum using a one and two tem- 
perature thermal plasma model (the MEKAL model in XSPEC) . 
A two temperature model gave a much better fit than a 
one temperature fit. The two temperature model (fcT=0.3, 
4.3keV) is better than the stratified accretion column model 
at the 94 percent confidence level. Again there is no evidence 
for a distinct soft X-ray component. The bolometric lumi- 
nosity (1 X 10^'^ ergs s~^ rfioo) is typical of systems in a low 
accretion state (=7.4x10^® ergs s~^for d=86 pc, Cropper 
1990). 
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0.016(3) 


0.33 


12203 


7392 
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9437 


MR Ser 


PG 1550+191 


571 


2003-01-22 


17 
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Table 1. The polars observed in our survey using XMM-Ncwtou wliicli were found to be in low accretion state or 'off' accretion state. The 
'Orbit' refers to the orbit number of XMA^-ATcwton and V was determined using the XMM-Ncwton Optical Monitor. We show the mean 
count rate derived using the EPIC pn with the exception of V834 Cen which was determined using the EPIC MOS camera. (For sources 
not clearly detected, we quote the background subtracted count rate for the expected source position). The first observation of V834 Cen 
was aifected by hardware problems. The RASS refers to the ROSAT all-sky survey (for AR UMa we used PSPC pointed observations, 
where it was not significantly detected, and for ST LMi where it was; Ramsay, Cropper & Mason 1995). The duration of the exposure 
in the EPIC MOS camera is shown in the last column. References for the V high state: (1) Imamura & Steinman-Cameron (1998), (2) 
Thomas & Reinsch (1996), (3) Singh et al (1995), (4) Burwitz et al (1997), (5) Schmidt et al (1999), (6) Imamura, Steinman-Cameron 
& Wolff (2000), (7) Reinsch et al (1994), (8) Thomas et al (1996), (9) Beuermann et al (1999), (10) Sekiguchi, Nakada & Bassett 1994 
(11) Cropper (1986), (12) Thomas et al (2000), (13) Motch et al (1996), (14) Burwitz et al (1996), (15) Liebert et al (1982). 
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Figure 2. The EPIC MOS2 spectra of V834 Cen. The best fit is 
shown as a solid line and is an absorbed two temperature thermal 
plasma model. 



3.2 UW Pic 

Reinsch et al (1994) present ROSAT survey observations of 
UW Pic and show that it has 2 maxima per cycle. Unfortu- 
nately their ephemeris is not precise enough to allow us to 
phase our data on their phase convention. However, we do 
find evidence for a variation in the X-ray light curve which 
is similar in shape to the ROSAT light curve, although the 
count rate is much lower (Figure 3). The signal to noise of 
the integrated spectrum is low. However, if we fit an ab- 
sorbed two temperature plasma model (fcT=0.5, 5keV) to 
the EPIC pn spectrum, we find an absorbed, bolometric lu- 
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Figure 3. The phased EPIC MOSl-l-2 data of UW Pic folded on 
the period of Reinsch et al (1994). 



minosity of 1 x 10^° ergs s~^for a distance of 300 pc (Rein- 
sch et al 1994). We find no evidence for a distinct soft X-ray 
component. 

3.3 V393 Pav 

Optical and X-ray observations reported by Thomas et al 
(1996) show V393 Pav to be a one-pole system with a bright 
phase lasting around 0.5 cycles. Based on the spectroscopic 
ephemeris the bright phase occurs at ~0.7-1.2. The accu- 
mulated error on the ephemeris up to the date of the XMM- 
Newton observations is 0.10 cycles. We show in Figure 4 the 
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Figure 4. The EPIC MOS (1+2) data of V393 Pav phased on 
the ephemeris of Thomas et al (1996). 



0.15-lO.OkeV light curve derived using XMM-Newton EPIC 
pn data. These data show a variable light curve with signifi- 
cant flux occurring at 4> ~0.0-0.4, although the bright phase 
in the second orbital cycle is more narrow than the first. 
Although the signal to noise of the spectra were low, we fit- 
ted them with an absorbed onc-tcmporaturc thermal plasma 
model (fcT=4.5kcV) and obtain a good fit and an unab- 
sorbcd bolomctric luminosity of 1.3x10"^" ergs s~^assuming 
a distance of 350 pc (Thomas et al 1996). There was no 
evidence for a distinct soft X-ray component. 



3.4 ST LMi 

ST LMi was not detected during the ROSAT all sky survey. 
Further, when it was the subject of a dedicated observation 
using ROSAT it was detected at the low count rate of 0.1 
ct/s (Ramsay, Cropper & Mason 1995). In contrast, the EX- 
OSAT observations of ST LMi show the system to be in a 
high accretion state and a bright X-ray target (ibid): it shows 
a bright phase lasting from 4> ~0. 7-1.1 on the ephemeris of 
Cropper (1986) which is coincident with the optical bright 
phase. 

We show the 0. 15-10. OkeV light curve made using EPIC 
pn data in Figure 5. It covers ~1.5 orbital cycles, with the 
expected maximum occurring near cf) ~0.0. On the first cycle 
no appreciable increase in flux is seen at this phase. How- 
ever, during the next cycle there is an increase in flux lasting 
A(f) ~0.2 cycles and centered on (f) ~0.0. There is no corre- 
sponding increase in the UVW2 count rate. 

We extracted a spectrum for the accretion event (which 
lasted ~1200 sec). We flnd it can be well modelled using 
an absorbed 2-tempcrature thermal plasma model (A:T=0.8, 
5keV) : no distinct soft X-ray component is required (Figure 
6). Using a distance of 128pc (Cropper 1990) we obtain an 
unabsorbcd, bolomctric luminosity for the accretion event 
of 1 X 10'*" ergs s^^. Extrapolating this best flt model to UV 



energies gives 4.3 x 10' 



-19 -1 

ergs s 



cm-^k in the UVW2 



filter. This is many orders of magnitude less than that ob- 
served (~ 1 X 10~^® ergs s~^ cm~^A). This suggests that 
the unheated white dwarf contributes a large fraction to the 
UV flux. 

If we assume that the X-ray event was due to accretion, 
we can estimate the mass of the accretion event through. 
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Figure 5. The phased EPIC pn data of ST LMi folded on the 
ephemeris of Cropper (1986). 



Figure 6. The EPIC pn spectrum of ST LMi fitted using an 
absorbed blackbody plus 2 temperature thermal plasma model. 



L = GMi M /Ri , where Mi and Ri are the mass and radius 
of the white dwarL Assuming Mi=O.7M0and Ri = 8.6 x 10* 
cm, M = 9.2 x 10^^ g s~^ or 1.1x10^*^ g for the duration of 
the event. This is remarkably similar to the mass of the ac- 
cretion event seen in UZ For (Still & Mukai 2001), although 
a corresponding increase in the UVWl band was seen in UZ 
For. This may imply a different temperature for the repro- 
cessed component in UZ For and ST LMi. 

3.5 RX J1313-32 

RX J1313-32 was discovered during the ROSAT all-sky 
survey when it was found to be bright reaching a maxi- 
mum count rate of 6.5 ct/s. When it was observed again 
using ROSAT in two pointed observations, it was found 
to be much fainter (0.04 ct/s in the PSPC, Thomas et al 
2000) . Our observations made using XMM-Newton show the 
system reaches a maximum count rate of 0.4 cts/s (EPIC 
MOSl+2 combined). Using the spectral flt described below 
we find using PIMMS (Mukai 1993) that we would get a count 
rate in ROSAT of 0.01 and 0.03 ct/s in the HRI and PSPC 
respectively. Although the system shows a low count rate 
it is still significant and it is likely that RX J1313-32 was 
observed in an low/intermediate accretion state (Figure 7). 

We show the spectrum taken using the EPIC pn detec- 
tor in Figure 8. There is prominent emission near the Fe Ka 
line. Closer inspection shows that this is due both to a line 
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Figure 7. The phased XMM-Newton data of RX J1313-32 folded 
on the ephemeris of Thomas et al (2000). Each bin size is 120s. 
The units for the OM data are 10~^® ergs s~^ cm~^A~'^. 
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Figure 9. The phased EPIC pn data of V1033 Cen folded on the 
ephemeris of CiesUnslii & Steiner (1997). 
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Figure 8. The integrated EPIC pn spectrum of RX J1313-32 
together with the best model fit. 



at 6.7keV and also the fluorescent line at 6.4keV. The spec- 
trum can be well fit (x^=l-04) using a stratified accretion 
column (Cropper et al 1999). Wc assumed a white dwarf 
mass of 0.7Moand 0.1 g s~^ cm~^. There was no evidence 
for a distinct soft X-ray component. The unabsorbed bolo- 

ergs s~ 

(Thomas et al 2000). 



metric luminosity is 5x10^° ergs s ^for a distance of 200pc 



3.6 V1033 Cen 

We show in Figure 9 the 0.15-lO.OkeV light curve of V1033 

Cen folded on the ephemeris of Cieslinski & Steiner (1997). 
It shows a variation in flux over the orbital period showing 
maximum near cf> ~ 0.6 and a minimum near <f) ~0.7-1.2. 
However, the accumulated error in the phasing is 0.37 cycles. 
Optical photometry shows a broad maximum around ~0.0 
(Buckley ct al 2000). Because of the phase error we cannot 
be certain whether the X-ray emission originates from this 
accretion pole. 

The spectrum has a low signal-to-noise ratio. How- 
ever, fitting a one-temperature thermal plasma model 
(fcT=3.5kcV) wc obtain a luminosity of 1.8x10^^ orgs s^^ 
dfoo for the integrated spectrum (the distance to V1033 Cen 
is currently unknown). There is no evidence for a distinct 
soft X-ray component. 
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Figure 10. The phased EPIC MOS(l-|-2) data of MR Ser folded 
on the ephemeris of Schwope et al (1993). 



3.7 MR Ser 

MR Ser was discovered using the Palomar Croon survey of 
blue stellar objects (Liebert et al 1982) and classified as a po- 
lar based on its optical spectrum and variable polarisation. 
Our observations of MR Ser show a clear variation in the X- 
ray flux (Figure 10). It shows a positive signal throughout 
the orbital cycle with a bright phase lasting ~0.4 cycles. It 
is reasonably similar to the light curve obtained using EX- 
OS AT on April 1 1984 (Angelini et al 1990). MR Ser was 
also found to be in a low accretion state during both the 
ROSAT all-sky survey and pointed observations. Wc know 
of no observations since then which have found MR Ser to 
be in a high state. We can obtain a good fit to its X-ray 
spectrum with an absorbed thermal plasma model with tem- 
perature ~1.4kcV. The unabsorbed bolometric luminosity is 
~ 3 X 10^'' ergs s"^ d' 
al 1993) gives - 6 x 10 



iQo, which for d=139 pc (Schwope et 
^® ergs s~^. 



4 SYSTEMS SHOWING NO X-RAY 
VARIATION 

4.1 CP Tuc 

Misaki et al (1996) found a new X-ray variable object which 
they classified as a polar. It showed a strongly modulated 
light curve which was most prominent at lower energies. Our 
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Table 2. Those systems which were detected in X-rays but found 
no X-ray modulation. We show the secondary mass estimated 
using the mass-period relationship, and the maximum X-ray lu- 
minosity assuming the mass-luminosity relationship of Malkov, 
Piskunov & Shipil'kina (1997) and a saturated X-ray luminosity 
of Lx /Ltoi = 1 X 10^'^. Wc also show the estimated distance and 
luminosity using the count rate to flux conversion described in 
the text. 

observations show that it was detected at the 6a level. How- 
ever, it does not show evidence for a significant modulation. 

4.2 MN Hya 

MN Hya was discovered during the ROSAT all-sky survey. 
Pointed ROSAT observations show a bright X-ray source 
with an absorption dip just before the ecUpsc of the sec- 
ondary star (Buckley et al 1998). The XMM-Newton ob- 
servations show a detection at the 5cr level. There was no 
significant variation in the light curve. 

4.3 The source of X-rays in these systems 

Is the source of X-rays in these systems due to accretion at 
a very low rate or is it due to emission from the secondary 
star? Late type stars show significant X-ray emission, gen- 
erally in the form of flaring activity. We now consider if the 
X-ray flux from these sources can be attributed to the sec- 
ondary star. 

Randich et al (1996) found that in late type systems 
Lx/Lboi does not exceed ~ 1 x 10~^. We show in Table 2 the 
secondary mass derived using the standard mass-period rela- 
tionship and the maximum X-ray luminosity using the mass- 
luminosity relationship of Malkov, Piskunov & Shipil'kina 
(1997) assuming a saturated X-ray flux. 

A late type secondary would show an X-ray spectrum 
dominated by a thermal plasma of temperature '^IkeV. To 
estimate the expected count rate, we assume this emission 
model and set the absorption to be 1 x 10^° cm~^. For the 
EPIC pn detector and a thin fllter, we find that a count rate 
of 0.01 ct/s would give a luminosity of ~ 2 x 10^^ ergs 

We find that in both CP Tuc and MN Hya the inferred 
luminosity is significantly greater than that expected for a 
late type main sequence star. We conclude that we have not 
detected the secondary star in either system and the flux is 
probably due to continuing low-level accretion. 



5 THE UV OBSERVATIONS 

Out of the 16 systems which were found to be in a low state 
using XMM-Newton, 10 were detected in both the UV fllters 
(UVWl has an effective wavelength of 2910 A and UVW2 
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Figure 11. We show the UVW1/UVW2 count rate ratio for dif- 
ferent white dwarf temperatures and absorption. We have used 
Hydrogen white dwarf model atmospheres (those of Koester) and 
the effective areas of each UV fllter. The count ratio for each low 
state system is marked at the left hand side. 

2120 A ). Another 3 systems (FH UMa, CV Hyi and RS Cae) 
were detected in the UVWl filter but not the UVW2 fllter. 
A further 3 systems, (AR UMa, V4738 Sgr and RX J0153- 
59) were detected in the UV but not in X-rays. UZ For had 
data taken in only the UVWl filter while V834 Cen had 
a technical problem with the UVW2 filter observation. We 
show in Table 3 the counts in each filter for the 10 systems: 
we have corrected the observed count rate for coincidence 
losses (Mason et al 2001). 

As part of the XMM-OM UV cahbration process, 8 
white dwarfs of various temperatures have been observed 
using the OM. These give conversion factors of count rate 
to fiux (1 ct/s in UVWl corresponds to a fiux of 4.4x10"^^ 
ergs s"^ cm-^A and 1 ct/s in UVW2 5.8x10"^^ ergs s"^ 
cm~^A- taken from the ESA XMM web site). These values 
give fiuxes to within 8 percent of their fluxes at these wave- 
lengths as observed using HST. We show the implied fluxes 
in the two filters in Table 3. 



5.1 Comparing UV count ratios with white dwarf 
models 

We show in Figure 11 the ratio of UVW1/UVW2 counts 
predicted by convolving white dwarf Hydrogen atmosphere 
models of different temperatures (those of Detlev Koester) 
with the effective area curves of the UV filters. We predict 
the ratios for a range in extinction (A^j^ — i5s_v'4.8 x lO'^^ 
cm~'^, Bohlin, Savage & Drake 1978). Polars typically have 
relatively low levels of extinction, so we take an 'upper-limit' 
of Eb-v = 0.1. We used the extinction law of Seaton (1979). 

The UVWl /UVW2 count ratios shown in Table 3 arc all 
greater than 5. Assuming low extinction this implies temper- 
atures lower than ~14000K. Those systems with the high- 
est ratios (CP Tuc and MN Hya) imply temperatures of 
~9000K for low extinction. These temperatures are signifl- 
cantly cooler than the temperatures inferred from UV spec- 
tra of polars in a high accretion state: typically an unheated 
white dwarf has T ~15000-25000K together with a hotter 
spot (Gansicke 1998). 
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5.2 Comparing low state and high state systems 

We now compare the UV ratios for systems in the low state 
with those polaxs which were observed in a high accretion 
state (these are taken from Ramsay & Cropper 2004). Table 
3 shows that the UVW1/UVW2 Count ratio for the low 
state systems is >5, with a mean of 7.4. In contrast the high 
state systems show a ratio <1.2 with a mean of 0.7. There 
is a clear difference between the high and low state systems: 
at face value this implies that in the high state the white 
dwarf is hotter compared with low state. 

However, in the high state it is common to have a soft 
X-ray component which is due to re-processing of hard X- 
rays from the photosphere of the white dwarf. This will 
make some contribution to the UV flux. To explore this fur- 
ther, we plot the ratios 0.15-0.5keV/UVWl (flux) against 
UVW1/UVW2 (flux) in Figure 12 for the systems shown in 
Table 3 and also those systems which were observed in a high 
accretion state. There is a clear separation of low and high 
state systems, with low state systems lying towards the left 
hand side, with high state systems in the right hand side. 

To determine if these ratios are consistent with that we 
expect, we estimate the UVW1/UVW2 and 0.15-0.5/UVWl 
flux ratios for both low and high state systems. For low 
state systems we take the best fit X-ray spectral model sim- 
ilar to that of ST LMi (§3.4): namely a two-temperature 
thermal plasma model. We also add a blackbody compo- 
nent of varying temperature which simulates the unhcated 
white dwarf (O.7M0and l.OMQand assume the Nauenberg 
1972 mass-radius relationship for white dwarfs) and place 
it at the distance of ST LMi (128pc Cropper 1990). (Since 
these estimates are already uncertain and we are using ra- 
tios, we use a blackbody rather than a white dwarf model 
atmosphere for the unheated white dwarf as is sufficient for 
our purposes). We assume the fraction of the white dwarf 
which is X-ray heated is much smaller than total surface 
area. For the high state systems we assume a model con- 
sisting of a soft X-ray blackbody component (of varying 
temperature) and a thermal bremsstrahlung component of 
temperature 30keV. Further, we assume an X-ray flux distri- 
bution consistent with the standard accretion model, ie the 
soft X-ray luminosity is approximately half that of the hard 
X-ray luminosity (cf Ramsay & Cropper 2003b). We also 
include a second blackbody which simulates the unheated 
white dwarf: we assume a temperature of 20000K (=1.7eV; 
appropriate for the white dwarf in a high accretion state, 
Gansicke 1998) and a 0.7Mq white dwarf. We scale the nor- 
malisation of this component so that the soft X-ray compo- 
nent due from re-processing covers a fraction of 10~* of the 
surface of the white dwarf. In both cases we fix the absorp- 
tion at 5x10^® cm~'^. For increasing levels of absorption the 
ratio 0.15-0.5keV/UVWl will decrease (we include a com- 
ponent for extinction in the UV based on the reddening law 
of Seaton 1979). 

We plot these tracks in Figure 12. For low state systems 
which have no reprocessed X-ray component we expect their 
ratio's to lie near the two left hand tracks: eg V1033 Cen, 
CP Tuc and MN Hya. Systems which have low temperature 
reprocessed components will lie to the lower right of these 
tracks: eg MR Ser, UW Pic, RX J1313-32, ST LMi and V393 
Pav. Despite the X-ray spectra of the low state systems not 
requiring a reprocessed X-ray component, our soft X-ray 
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Figure 12. We show the observed 0.15-0.5keV/UVWl and 
UVW1/UVW2 flux ratios for those systems in a low accretion 

state and a high accretion state (taken from Ramsay & Cropper 
2004). For high state systems we show tiie expected ratio assum- 
ing a polar with mass O.TMq (and using the Nauenberg 1972 
relationship) and temijcrature 1.7eV (=20000K) and which has 
a spectrum consistent with that predicted by the standard ac- 
cretion shock model (we show the ratio for various temperatures 
for the reprocessed soft X-ray component). We also show the ex- 
pected ratio for low state systems assuming a spectrum like that 
of ST LMi: we also add a blackbody of various temperature to 
account for the unheated white dwarf. We assume the interstellar 
absorption to be 5x10^" cm~^. 

and UV data imply that even in the low state, a reprocessed 
component exists: it is not detected in X-rays because the 
temperature is low enough to have moved to energies lower 
than that of soft X-rays. In this case, the UVW1/UVW2 
ratio is insensitive to the temperature of the unheated white 
dwarf until temperatures of less than ~5eV are reached. 



6 THE FREQUENCY OF SYSTEMS IN A LOW 
ACCRETION STATE 

In our HMM-l^ewton X-ray survey, we observed a sample 
of 37 systems, which comprises a significant portion of the 
currently known polars. The orbital periods of this sample 
range from 77 min to 270 min. Sixteen systems (43 percent) 
of the systems were found in a state of low or much reduced 
accretion rate. 

Mass transfer in magnetic cataclysmic variables, includ- 
ing polars, is believed to be governed by the orbital evolution 
of the binary The frequency of high/low state transitions 
and the durations of the state duty cycles determine the ob- 
served proportion of high or low-state systems. Wo ask: are 
the high/low state transition and duty cycle of polars gov- 
erned by processes dependent on the orbital period? If so, 
association of the high/low states with the orbital periods 
will be reflected in our XMM-JVewton polar survey sample. 
To verify this we carry out several statistical tests, (i) to 
examine the bias of the high and low-state systems toward 
long or short orbital periods, (ii) to verify whether or not 
there is an alternating sequence of high/low state transition 
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Source 




UVWl 




UVW2 


0.15-0.5keV 


UVW1/UVW2 




Count 


Flux 


Count 


Flux 


Flux 


(Counts) 




Rate 


xlO-i'' 


Rate 




xlO-i-'^ 






(Ct/s) 


ergs s^-"^ cm~^A 


(Ct/s) 


ergs s^-"^ cm^'^A 


ergs s~^ cm~^ 





CP Tuc 


0.60(3) 


2.6(1) 


0.06(1) 


3.5(2) 


1.7 


10.0 


AR UMa 


10.08(9) 


44.4(5) 


1.61(5) 


93(2) 




6.3 


UW Pic 


1.62(7) 


7.1(1) 


0.25(2) 


14.5(5) 


10.0 


6.5 


V393 Pav 


0.53(3) 


2.3(1) 


0.07(1) 


4.1(6) 


7.0 


7.6 


V4738 Sgr 


0.15(3) 


0.7(1) 


0.03(1) 


1.7(6) 




5.0 


RX J0153-59 


0.42(2) 


1.9(1) 


0.08(1) 


4.7(7) 




5.3 


MN Hya 


0.36(3) 


1.6(1) 


0.03(1) 


1.7(2) 


3.5 


12.0 


ST LMi 


1.78(3) 


7.8(1) 


0.25(1) 


1.5(6) 


35 


7.1 


RX J1313-32 


5.81(8) 


26.1(3) 


1.10(3) 


64(1) 


55 


5.3 


V1033 Cen 


0.61(3) 


2.7(1) 


0.07(1) 


4.1(6) 


1.7 


8.7 


MR Ser 


3.11(4) 


14.0(2) 


0.40(2) 


23(1) 


16 


7.8 



Table 3. For those systems which were detected in both the UVWl (peak response 2910 A) and UVW2 (21210 A) filters we show their 
count rate in each filter and the corresponding flux. The number in parenthesis is the error (standard deviation) on the last figure. We 
also show the observed fiux in the 0.15-0.5keV energy band. 



(including periodic transition) which can introduce clump- 
ing of particle states at certain orbital periods, and (iii) to 
determine the relative duration of high/low-state duty cy- 
cles. 

We note that polars have previously been surveyed us- 
ing ROSAT. We therefore also re-analyse the ROSAT sample 
and compare it with our XMM-Newton sample. In defining a 
suitable sample from the ROSAT data, we consider all polars 
which were not discovered using the ROSAT all-sky survey. 
We also exclude the very faint systems which may not be de- 
tected in the relatively short exposure of the ROSAT all-sky 
survey (eg V1500 Cyg, EU Cnc which are relatively distant 
systems). Wc performed a similar study to that done for the 
XMM-Ncwton sample - namely wo compared the count rate 
for each source determined during the ROSAT all-sky survey 
with other X-ray observations of the same source. Out of a 
total of 28 systems, 16 (57 percent) were in a low accretion 
state during the ROSAT all-sky observation. 

In the statistical analysis of the two samples, we assume 
that the duration of the high state and low states (^ weeks, 
months or longer) are much longer than exposure time of 
individual systems (typically <10 ksec) during the survey. 
Under this assumption, the observation of each system can 
be treated as a snapshot. 




Orbital Period (mins) 

Figure 13. The period distribution of polars in high and low 

accretion state (assigned with values 1 and respectively) during 
the XMM-Ncwton survey (top panel) and the ROSAT all-sky- 
survey. We also mark the periods of 120 and 180 min (dotted 
vertical lines), the approximate boundaries of the period gap, the 
range of periods where magnetic cataclysmic variables have been 
claimed to be deficient. 



6.1 High/low state duty cycle and orbital-period 

As we do not have information about the functional forms 
of the orbital period distributions of the high-state and the 
low-state systems, conventional parametric analyses such as 
those involving the means and the variances of the two dis- 
tributions are not directly applicable. To overcome this diffi- 
culty, wo consider non-paramotric tests instead. These tests 
can be applied even for situations in which the numerical 
scales of the measurements are arbitrary, thus allowing us 
to examine how the high/low state transition probability is 
related to the orbital period in an objective manner, without 
pre-assumed models for the association. 

We consider a Wilcoxon-Mann- Whitney rank sum U- 
test (Wilcoxon 1945, van der Waerden 1969) to search for 
biases of high and low-state systems towards long/short or- 



bital periods and a general tost of randomness (see eg John- 
son 1994) to examine the possibility of orbital clumping of 
high-state and low-state systems. Both tests require the sys- 
tems to be ranked, and we choose a simple criterion: rank- 
ing according to the orbital period. We also adopt the two 
following conventions in our analysis. (1) The system with 
the shortest period has the 1st rank, and a higher rank for 
longer-period systems. (2) Each low-state system is assigned 
with a state identity parameter of value L, each high-state 
system with a state identity parameter of value H. The ranks 
of the systems in the XMM-JVewton and ROSAT surveys are 
shown in Table 4, while the period distribution of the sys- 
tems are shown in Figure 13. 
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XMM-Newton 



ROSAT 



proximately by a normal distribution, with the mean and 
variance given by 
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77.8 
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J0132-65 


80. 
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FH UMa 


78.0 


L 


J2022-39 


81.02 


H 


EF Eri 


79.7 


H 


EV UMa 


85.8 


H 


J0425-57 


79.9 


H 


J1015+09 


89.4 


L 


CP Tuc 


80.0 


L 


J0154.0-5947 


89.8 


L 


DP Leo 


80.15 


L 


FH UMa 


100.4 


H 


VV Pup 


89.04 


L 


CP Tuc 


101.5 


H 


V834 Cen 


89.8 


H 


DP Leo 


104.6 


H 


EP Dra 


90.08 


H 


V347 Pav 


108.6 


L 


CE Gru 


90.14 


H 


J1149+28 


112.97 


H 


V2301 Oph 


95.0 


L 


J0453-42 


113.47 


L 


MR Ser 


98.8 


L 


J1957-57 


113.64 


H 


BL Hyi 


101.5 


L 


V834 Cen 


113.9 


L 


ST LMi 


104.6 


H 


EP Dra 


114.6 


L 


WW Hor 


107.0 


H 


J1002-19 


114.84 


H 


AN UMa 


108.6 


H 


CE Gru 


115.45 


H 


EK UMa 


110.9 


H 


J2115-58 


115.92 


L 


AR UMa 


113.6 


L 


MR Ser 


126.5 


L 


UZ For 


113.9 


L 


ST LMi 


156. 


L 


SDSS1324 


114.84 


H 


AN UMa 


160. 


L 


V349 Pav 


115.45 


H 


EK UMa 


167.7 


L 


WX LMi 


115.49 


H 


WW Hor 


185.6 


L 


AM Her 


115.92 


L 


AR UMa 


199.3 


H 


BY Cam 


125.4 


H 


EU One 


222.5 


H 


QQ Vul 


126.5 


L 


UZ For 


228.4 


L 


VY For 


128.7 


H 


J1846+55 


244. 


L 


HS0922+13 


133.0 


L 


J0531-46 


263.4 


L 


SDSS1553 


160.0 


H 


V2009-65 


270. 


L 


V895 Cen 


171.3 


H 


J0501-03 








189.4 


L 


J1141-6410 








199.3 


H 


BY Cam 








201.0 


H 


VI 500 Cyg 








203.4 


L 


MN Hya 








208.0 


H 


J1007-20 








228.4 


H 


VY For 








251.4 


L 


J1313-32 








270.0 


H 


V895 Cen 









Table 4. The list of sources together with their orbital period 
and flagged for whether they were in high or low accretion states 
at the epoch of our XMM-Newton survey (left) or the ROSAT 
all-sky survey (right). For reasons of space we have omitted the 
'RE' and 'RX' prefix for sources discovered using ROSAT. 



m(C/i) = KU2) 



nin2 



nin2(ni + n2 + 1) 
12 



The associated standardised random variable is defined as 
Ui - njUi) 

Thus, the statistics that characterise the XMM-Newton 
sample are Ui = 132, U2 = 204 (cf Table 4), n{Ui) = 168 
and a'^{Ui) = 1064, implying a standardised random vari- 
able «xmm = —1.10. 

In the ROSAT sample, we have ni = 16 and 7x2 = 11, 
so the U statistic can again be approximated by a normal 
distribution. The rank sums of the two types of systems 
are W\ = 279 and W2 = 127; the relevant statistics are 
Ui = 143, U2 = 49, ^{Ui) = 96 and cr'^iUi) = 464. The cor- 
responding standardised random variable is Zrosat = +2.18. 

For 99 and 95 percent confidence interval, —2.575 < 
z < 2.575 and —1.96 < z < 1.96 respectively. The z value 
for the XMM-Newton sample is well within the 95 percent 
confidence interval, while the ROSAT sample is within the 
99 percent confidence interval. Therefore, we cannot reject 
the null hypothesis on a high significance level to favour the 
alternative hypothesis. 

We note that similar results can be obtained by the 
Kolmogorov-Smirnov test, which may be inefficient to dis- 
criminate the difference in the period distribution of the 
high-state and low-state systems especially when the dif- 
ference occurs near the two ends of the period range. The 
D-statistic that we obtain for the XMM-Newton sample is 
0.277, corresponding to a probability of 0.49 for the null 
hypothesis. The ROSAT sample has a larger D-statistics, 
0.479, and a null-hypothesis probability of 0.09. 

As a summary, we conclude that the XMM-Newton 
sample and the ROSAT sample do not show significant dif- 
ference in the (global) distributions of orbital periods among 
the high-state and low-state systems between 77 and 270 
min. That is, there is no firm bias of high-state or low-state 
systems toward short orbital periods and towards long or- 
bital periods. 



6.1.1 Bias in the orbital periods 

First, we verify if there is an excess of either high-state or 
low-state systems among the long-period polars and among 
the short-period polars. The null hypothesis is: there is no 
bias of high-state and low-state systems towards long orbital 
periods or short orbital periods. 

In the XMM-Newton sample, the number of low-state 
systems is rii = 16 and the number of high-state systems 
is 712 = 21. The corresponding rank sums of the two types 
of systems are W\ = 268 and W2 = 435. In the Wilcoxon- 
Mann- Whitney test, the [/-statistic is given by 

U, = H^,-M!|±l) (, = 1,2) 

As the numbers of two types of systems are sufficiently large 
(ni > 10 and n2 > 10), the ZJ-statistics are described ap- 



6.1.2 Period clumping of high and low-state systems 

We now test whether the systems in the two accretion states 
are in a random order when ranked according to their or- 
bital periods. This test will reveal the sequence of alternating 
high/low states in the population if the alternation depends 
strongly and explicitly on the orbital period, but it does 
not exclude such transitional behaviour for each individual 
system. 

The null hypothesis is: the period ordering of the high- 
state and low-state systems are random. Using our adopted 
convention, we construct a sequence for the XMM-Newton 
sample: LLHHL... (cf Table 4). We then define a w-statistic, 
which is the number of runs (sec .Johnson 1994) in the se- 
quence, and for the XMM-JVewton sample it is 20. For sam- 
ples with m > 10 and n2 > 10, the w-statistics approxi- 
mately have a normal distribution. The mean and variance 
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of the distribution given by 
2nin2 



2nin2(2nin2 — ni 



■ n2 



(ni + n2)2(ni + n2 - 1) 
the standardised random variable of the w-statistic is 

u — i-i{u) 
a{u) 

The statistics of the XMM-Newton sample are therefore 
u = 20, = 19.16 and o-^(w) = 8.66. It follows that 

the standardised random variable ^xmm = 0.29. 

Similarly, we can order the ROSAT systems in a similar 
manner. This gives u = 12, fj.{u) = 13.71, cr^{u) — 6.46 and 
2:rosat = —0.67. As Zxmm and «rosat IS Well within the interval 
— 1.00 < z < 1.00, the null hypothesis is consistent with the 
data of both survey, implying that the period ordering of 
the systems are probably random. We note that visual in- 
spection may suggest a deficiency of high-state systems with 
periods between 120 and 180 min (the controversial period 
gap of polars) in the ROSAT sample (Figure 13). However, 
this apparent deficiency of high-state systems is not seen 
in the XMM-Newton sample. The randomness test that we 
have performed indicates that any dumpiness in the two 
sequences can be explained by statistical fluctuation. Com- 
bining the result of this test and that in §6.1.1, we conclude 
that the XMM-JVewton and ROSAT data do not show evi- 
dence that the high-state and low-state transition and duty 
cycle of polars are associated with certain preferred orbital 
periods. 



6.2 Probability of being in the low state 

Inspection of the proportion of low-state sources suggests 
that the fractional duty cycle of low accretion state is 0.5. 
However, what is the uncertainty in this fraction? We show 
in the Appendix how we calculate this uncertainty in a for- 
mal Bayesian approach. 

In the XMM-Newton sample N = 37 and k = 16, and 
in the ROSAT sample iV = 28 and k = 16, where iV is the 
number of systems in the sample and k is the number of 
systems in a low accretion state. Figure 14 shows the cu- 
mulative probability distribution for low-state duty cycle x, 
which is P{[x, 1] \k, N), and its derivative for the two surveys. 
The estimate x is 0.442 and the 90% confidence interval is 
(0.258, 0.734) for the XMM-Newton sample. The estimate x 
is 0.585 and the 90% confidence interval is (0.342, 0.863) for 
the ROSAT sample, consistent with the result obtained from 
the XMM-Newton sample. Thus, we may conclude that the 
probability of a system in a low state is ~50 percent, with an 
uncertainty corresponding to the two samples given above. 



6.3 Implications for the space density of polars 

Three new polars were discovered using EXOSAT in a cov- 
erage of 5 percent of the sky. Based on this finding Hameury, 
King & Lasota (1990) predicted that around 60 new polars 
would be discovered using ROSAT. Since ROSAT was more 
sensitive than EXOSAT), this number would be expected 
to be greater. In contrast, ~30 polars were discovered using 
ROSAT. However, the Poisson uncertainty of the prediction 
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Figure 14. The cumulative probability distribution P{x; k) = 
P([a;, 1] I k, N) (solid line) for the low-state duty cycle x and its 

derivative dP/dx (dotted line) for the ROSAT a.nd XMM-Newton 
surveys (top and bottom panels respectively). The best estimate 
of the low-state duty cycle, x, is given by the peak of dP/dx, 
and the bounds of the 90% confidence interval are located by 
P{x; k) = 0.95 and 0.05 respectively. 



is ±35. Taking into account of statistical fluctuations, we 
conclude that the detection of approximately 30 systems by 
ROSAT is consistent with the results of the EXOSAT sur- 



vey. 



10 



The estimated space density of CVs is ~ (1 — 3) x 
pc"^ (Patterson 1998, Schwope et al 2002a). The polar 



to disk CV ratio, inferred from the Einstein galactic plane 
survey, is about 1/3 (Hertz et al 1990). Thus, the space den- 
sity of polars is probably slightly less than 10~^ pc~^. 

Because of inhomogeneities in the local interstellar 
medium, soft X-rays sources have very different discovery 
probabilities in different parts of the sky. Soft X-ray sources 
such as polars and single white dwarfs, could be preferen- 
tially found in certain parts of the sky (eg Warwick et al 
1993). Thus, the survey carried out by the EXOSAT LE 
inevitably suffered from this bias. 

In addition to the bias duo to absorption by an inlio- 
mogeneous local interstellar medium, there are several other 
limiting factors that could also introduce uncertainties in 
the estimates of the space density of galactic polars. Among 
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them is the lack of accurate distance estimates for many 
systems (many polars have only lower limits), and the fact 
that those polars which accrete at an extremely low rate (cf 
Schwope et al 2002b) are unlikely to be discovered in X-ray 
surveys. Another important limiting factor is the uncertain- 
ties in the proportion of high-state and low-state systems 
in the population. In deriving the space density for polars, 
many workers (eg, Thomas & Beuermann 1998) first as- 
sumed that the sample is complete at high galactic latitudes, 
then corrected for parts of the sky (such as the North-Polar 
Spur) which have higher interstellar absorption and finally 
rescaled according to the fraction of systems in the high- 
accretion states. This survey has eliminated one uncertainty 
by showing that about a half of the polar population is in a 
low accretion state at a given time. 



7 CONCLUSIONS 

In our survey of 37 polars using XMM-Newton we found 

that 16 were a low accretion state. Of these, 6 were not de- 
tected in X-rays. The XMM-Newton data show no evidence 
that the low state systems are biased towards certain or- 
bital periods. The period ordering of the system is consistent 
with random. Similar results were found when re-examining 
a sample of the ROSAT all-sky survey data, in which 16 out 
of 28 systems were in low accretion states. 

Of those low state systems which were detected at a 
significant level, 8 showed a significant variation in X-rays. 
We find that in epochs with low accretion states, accretion 
is still occurring, often sporadically at the same accretion 
pole as when the system is in a high accretion state. Their 
spectra can be modelled using a two-temperature thermal 
plasma model of temperature several keV. The unabsorbed 
bolometric luminosity is typically ^ 10^° ergs s~^: this is 
two orders of magnitude less than the luminosity observed in 
the high accretion state (Ramsay & Cropper 2003b). There 
is no evidence for a distinct soft X-ray luminosity as often 
seen in the high accretion state. However, the UV data sug- 
gests that a reprocessed component is still present in the 
low accretion state, but is cooler than in the high state and 
therefore shifted into the UV or far-UV. 
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APPENDIX: Determining the probability of being 
in the low state 

The probability of m systems in a low state in a sample of 
N system is given by a binomial distribution 



b{m,N;x) = 



m\{N--^ 



where X is the low-state duty cycle of the systems. If the sys- 
tems are not observed simultaneous, it introduces counting 
fluctuations, which could be approximated by a Poisson pro- 
cess. The probability of k systems in a sample of A'^ systems 
are observed in the low state during a survey is therefore 



a 



E 



X / dxx"'{i-xy 



k\ m\{N-my. 



Bb{m + 1,N -m+1) - Ba{m + \,N -m + l) 



where B),{m -|- 1, A'" — m -|- 1) and Ba{m + 1, N — m + 1) are 
incomplete Beta functions. The denominator 7(0, 1; k, N) is 
given by 



7(0,1; A;, AT) = ^ 



= E 



k\ m\{N-my. 

m=l 

X / da; a;'"(l - x)^-"" 
Jo 

m T{m + l)r{N -m+l) 



k\ m\{N-my 



r{N + 2) 



+ 1) E 



{N+1) ^ k\ 

m=l 

where T{m + 1), T{N -m + l) and r{N + 2) are Gamma 
functions. 



P{k,N\x) = ^ 



m— 1 

JV 



A;! 



= E 



6(m, N; x) 



k\ m\{N-my. 



X (1 — a;) 



From the Bayesian Theorem, the probability of a system 
being in a low state given that k systems are observed in a 

low state is 

f'dx P(k,N I x)f(x) 

P(a<x<b\k,N) = — L_LL^ 

J°dxP{k,N I x)f{x) 

where f{x) is a priori probability density distribution of x. 
For a uniform density distribution f{x), we have 



P{a<x<b\k,N) = 



I{a,b-k,N) 



7(0,l;fc,iV) ' 

and X can be interpreted as the low-state duty cycle. The 
numerator 7(a, 6; k, N) is given by 



I{a,b;k,N) 



N 
m=l 



N] 



k\ m\{N-my 
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